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Abstract. The recent distance determinations by HIP- 
PARCOS have allowed to locate very precisely a reason- 
able set of subdwarfs in the HR diagram. A detailed com- 
parison with evolutionary models is now possible and it 
has been recently claimed that the observed lower part 
of the main-sequence is cooler than the computed one. In 
this paper we show that microscopic diffusion could be 
the physical explanation of such a systematic shift. We 
emphazise the fact that, as diffusion alters surface abun- 
dances, the initial metal content of the stellar material is 
higher than the present one. A "calibration" procedure is 
necessary when comparing observed and theoretical HR 
diagram, to take into account this effect. 
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K*" ■ 1. Introduction 



Microscopic diffusion is a basic phenomenon, but, in the 
physical conditions of stellar interiors the small diffusion 
velocities imply long time scales to obtain significant mo- 
difications on the global parameters of real stars. To be ef- 
ficient, the medium has to be quiet enough, so that large 
scale motion cannot prevent the settling. Acting essen- 
tially in radiative zones, it is now been advocated to ex- 
plain a variety of observed facts, as lithium abunda nces i n 
young and old population (Vauclair & Charbonnel |l99S|) , 
abundances anomalies i n diffe rent classes of young stars 
(Hui Bon Hao & Alecian 1998 and references herein) , glo- 
bular cluster evolution and age determination (Chaboyer 
et al. 1992 ). In the Sun, the insertion of microscopic diffu- 
sion in the modeling improves significantly the agreement 
between the theoretical and t he "se ismic" model (se e e.g . 
Christensen-Dalsgaard et al. [1996 ; Richard et al. 1996 ; 
Morel et al. 1997| ; Brun et al. |1998|) . In the present state 
of art the computed sound speed, beneath the solar con- 
vection zone, agrees within a rms discrepancy better than 



0.2% (Gough et al. 1996); the predicted values of the ra- 
dius at the bottom of the convection zone and the helium 
abundance at the solar surface agree within the error bar 
with their values inferred from helioseismology (see e.g. 
Basu[l997D. 



Subdwarfs represent a clue in the general framework 
of stellar evolution in particular to estimate the globular 
clusters ages. A precise knowledge of their internal struc- 
ture and of the stage of evolution is then badly needed. 

As in such old objects the evolutionary time scale be- 
comes of the same order as the diffusion one, significant 
effects due to this physical process are expected on their 
global properties. 

Several papers have already mentioned the infiuence 
of mic rosco pic diffusion in old ob jects (see e.g. Ch aboye r 



et al. 1992; MaziteUi et al. 1995; Castellani et al. 1997), 



but focussed generally on the consequences on the age of 
the oldest globular clusters. 

Up to now, the poor quality of the observed parameters 
of population II objects did not allow to refine their the- 
ory. The recent distance determinations of subdwarfs by 
HIPPARCOS (Perryman et al. |l997a| ) has revolutionized 
the field. Combining these results with large improvement 
of effective temperature and bolometric corrections for low 
metal atmospheres (Alonso et al. 1996), a reasonable set 
of subdwarfs have been located very precisely in the HR 
diagram. A detailed comparison with evolutionary mod- 
els is now possible and requires the same precision in the 
description of the physics of their interiors. 

It has been recently emphasized that difficulties arise 
when trying to fit observed and computed main-sequences, 
as models look hotter than the real objects (Baglin 1997; 
Lebreton et al. 1997; Cayrel et al. 1997; Lebreton et al. 



1998) 



Several effects can been advocated to explain this dis- 
crepancy, as for instance, the treatment of the supera- 
diabatic outer layers or observational bias due to NLTE 
effects, in the effective temperature scale and in the abun- 
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dance determinations (Thevenin & Idiart 1999). 

We propose in this paper to document the role of mi- 
croscopic diffusion in these stages of evolution, and to 
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quantify its influence on the global parameters of these 
stars and on their position in the HR diagram. 

It is well known that, on the main-sequence, variations 
of helium Y and of the metal content Z act in opposite 
directions : while decreasing Y reduces the effective tem- 
perature Teff and the luminosity L, a decrease of the metal 
content increases both L and Toff. Microscopic diffusion, 
which acts on both Y and Z, creates a stratification of 
chemicals in radiative zones. The global effects on the ob- 
servables i.e. the modification of the HR diagram position, 
but also the changes of the surface metal abundance, are 
quite subtle; they have to be taken into account altogether 
to perform a precise comparison with observations. 

The paper is organized as follows : the set of stellar 
models is listed in Sect, g, whereas the input physics is 
described in Sect. |^. The effects of microscopic diffusion 
on the evolution of subdwarfs are presented in Sect. ^, the 
" calibration" procedure needed when comparing theoreti- 
cal predictions and observations is developped in Sect, ra , 
whereas the influence on the HR diagram and on the shape 
of the isochrones as a function of their present observed 
metal content is described in Sect. 0. Future prospects and 
possible tests of this hypothesis are suggested in Sect. |^. 

2. The stellar models 

As subdwarfs belong to population H, their initial helium 
content is close to the primeval value. Then stellar models 
depend only on age and metal content, assuming that the 
physical description is completely settled. 

In old metal poor stars, a elements enrichment seems 



Table 1. [Fe/H] and initial abundances, per mass unit, of 
hydrogen X and heavy elements Z. 



to be always present (see e.g. Wheeler et al. 1989), though 
large variations exist from one star to the other. We 



have used in all our models the mixture of AUard ( [1996D 
with the a-enrichmcnt relative to the Sun, [a/Fe] ~ 0.3. 
Four initial metal contents, represented by the stan- 
dard [Fe/H] observable, have been considered : [Fe/H] r^ 
-0.75, -0.94, -1.24 and - 1.72. The initial total amount 



^He, per mass is fixed to its primor 
(Izotov et al. |l997D 



of helium, i.e. He 

dial value assumed to be Y = 0.24 

With respectively Zpc, X and Z as the iron, hydrogen and 

heavy elements mass ratio, we use the relation : 



[Fe/H] = log(ZFc/Z) + log(Z/X) - log(ZFe/X)0 



(1) 



to relate the observed metallicity and metal content. Here 
log(ZFo/X)o = —2.753 for the solar mixture of Grevesse 
& Noels, (|1993D, and log(ZFc/Z) = -1.371 for the AUard's 



a-enriched mixture. 

The initial mass ratio X and Z are derived from the 
initial values of Y and [Fe/H] as given in Table U 

To describe the lower part of the main-sequence we 
present here models of masses varying from 0.6 to O.85M0 
evolving up to 10 Gy. Below 0.6Mq, evolution is extremely 
slow, and above 0.85Mq subdwarfs have left the main- 
sequence. As during the rapid pre main-sequence phase, 
the stellar models are almost fully mixed by convection. 



Labels [Fe/H] X 



07 -0.75 0.7544 0.0056 

09 -0.94 0.7564 0.0036 

12 -1.24 0.7582 0.0018 

17 -1.72 0.7594 0.0006 



the microscopic diffusion is inefficient, and starting the 
computations at the homogeneous main-sequence stage is 
sufficient for our purpose. 

3. Physics of stellar models 

Nuclear network and chemical mixture. The general nu- 
clear network we used contains the following species : "'^H, 
^He, -^He, ^^C, ^^G, "N, ^^n, isq, "O, and Ex. Ex is 
a fictitious mean non-CNO heavy element with atomic 
mass 28 and charge 13, which complements the mixture. 
With respect to time, due to the diffusion processes, the 
abundances of heavy elements are enhanced toward the 
center; Ex mimicks that enhancement for the non CNO 
metals which contribute to changes of Z, then to opacity 
variations. We have taken into account the most impor- 
tant nuclear reactions of PP-f CNO cycles (Clayton, 1968 ) 
with the species ^H, ^Li, ^Be at equilibrium. The relevant 
nuclear reaction rates are taken from Caughlan and Fowler 
( 1988 ); weak screening is assumed. 

For helium the initial isotopic ratio is fixed at 
^Re/^Re= 4.19 lO"'' (Gautier & Morel |1997D . In the 
mixture of heavy elements, the ratios between the CNO 
species are set to their a-enriched AUard (1996) val- 
ues (in number) C : 0.147909 , N : 0.038904 and O : 
0.616594, with the isotopic ratios : ^^C/^'^C^ 1.11 10"^, 
i5n/14n= 42510-3, 1^0/1*^0= 3.8110-'' (Anders & 



Grevesse, 1989) 



Diffusion. Different processes are participating to ele- 
ment separation or mixing. Whereas gravitational set- 
tling creates stratification, hydrodynamical instabilities 
can, at least in some phases of evolution, generate macro- 
scopic motions which tend to reduce the chemical inhomo- 
geneities. Presently, no general description exists of these 
processes, which could be easily incorporated in stellar 
evolution calculations. As the purpose of this paper is to 
evaluate the influence of elements segregation in compar- 
ing theory and observations of subdwarfs, we did not take 
into account any of these processes. Let us note that ro- 
tation, which has been identified as one major cause of 
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Fig. 1. Profiles of abundances, in mass, of hydrogen X 
(full), helium Y (dash), heavy elements Z (dot-dash) and 
mean molecular weight [i (dot) for two O.85M0 models 
evolved up to 10 Gy, respectively with (heavy) and with- 
out (thin) microscopic diffusion. 



partial mixing is almost absent in these old low mass ob- 
jects. 

Microscopic diffusion is described using the formal- 
ism of Michaud & Proffitt (1993) valid for main-sequence 
stars. The radiative forces are not taken into account. For 
Y the mass conservation equation gives : Vy = —XVx/Y 
where the original equation has been slightly modified to 
take into account the heavy elements; Vx and Vy are 
respectively the diffusion velocities of hydrogen and he- 
lium. This description is parameter free, depending only 
on known coefficients. A more complete discussion of the 
metho d used to treat diffusion is given in Morel et al. 
(119971) . 



Equation of state and Opacities. We have used the EFF 



equation of state (Eggleton et al. 1973 ) sufficient for our 
purpose in this mass range. 

We have chosen the Livermore Library (Iglesias & 
Roger s |1996 ) with the a-enriched mixture of AUard 
( |l996| ), complemented at low temperature opacities by the 
Kurucz's (1998) a-enriched tables. Unfortunately there is 
no available common mixture in these libraries; we have 
retained the closest ones with a rough smoothing in the 
extreme low temperature atmospheric layers. 

The opacities and equation of state are functions of 
the heavy elements content Z, through the number of free 
electrons and the abundances of efficient absorbers which 
do not necessarily belong to the nuclear network e.g. '^^Fe. 
Due to diffusion and nuclear reactions, Z changes as the 
star evolves, as well as the ratios between the abundances 
of chemicals. In the calculation of opacities and equation 
of state, Z is separated in two parts : the first one consists 
of the chemicals heavier than helium which, belonging to 
the CNO nuclear network, are both diffused and nuclearly 



processed; Ex, the second part, is only diffused. Hence in 
the estimate of Z, the changes of CNO abundances caused 
by diffusion, nuclear processes and the effects of the gravi- 
tational settling of the heaviest non-CNO species are taken 
into account. 



Convection. In the convection zones the temperature gra- 
dient is computed according to the standard mixing-length 
theory, with the mixing- length defined as I = XHp, where 
Hp is the classical pressure scale height. The constancy 
of A has been demonstrated for population I m odels close 
to the main-sequence (see e.g. Fernandes et al. 1997), and 
calibrations of A from 2D simulations (Ludwig et al. 1998 ) 
lead to approximately the same result; so, we assume that 
A is constant and equal to the value obtained for a so- 
lar model using the same physics, namely A = 1.7. This 
hypothesis has no influence on the results of this paper 
devoted to an estimate of the differential effect of micro- 
scopic diffusion. 

In the models with diffusion the convection zones are 
mixed via a strong turbulent diffusion coefficient, which 
produces an homogeneous composition. 



Atmosphere. An atmosphere is restored using the Hopf's 



T{t) law : T^ = IT^^[t + (j(t)] (Mihalas |1978| ) with r as 



the Rosseland optical depth. The connection with the en- 
velope is made at the optical depth Tb — 10, where the dif- 
fusion approximation for radiative transfer becomes valid. 
In the convective part of the atmosphere, a numerical trick 
(Henyey et al. 1965) is employed in connection with the 
purely radiative Hopf's law in order to ensure the con- 
tinuity of gradients at the limit between the atmosphere 
and the envelope. At each time step, the radius i?* of the 
model is taken at the optical depth r* ~ 0.6454 where 
T{Ti,) — Tcff ; the mass of the star M^, is the mass inside 
the sphere of radius i?-* . The external boundary is fixed at 
the optical depth Text = 10~^, where a boundary condition 
on the pressure is expressed as ; P(roxt) — 9 / i^Tcxti with 
g as the gravity and k as the Rosseland mean opacity. 

Numerics. The mode ls ha ve been computed using the 
CESAM code (Morel |l997|) . Typically each evolutionary 
track needs of the order of 85 models of about 500 mass 
shells. The accuracy of the numerical scheme is one for 
the time and three for the space. 

4. EfTects of microscopic diffusion on the 
structure of main-sequence models 

Due to the gravitational settling, all species but ^H sink; 
therefore as time goes on, the amount of hydrogen at the 
surface is enhanced while the amount of helium and heavy 
elements decreases. Then a first effect of microscopic dif- 
fusion is the decrease of the surface value of [Fe/H] with 
respect to time. The variation of Z modifies the opacity 
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Fig. 2. The same, see Fig. ^ for normalized pressure 
(full), temperature (dash) and luminosity (dot-dash). The 
factors of normalization are respectively 4.710^^ dyne 
cm~^ for the pressure, 1.710^ K for the temperature and 



3.6 10^3 erg cm^^ g-i 



for the luminosity. 



and then changes slightly the temperature gradient in the 
radiative zones, as well as the position of their boundaries. 
In the mean time, the sinking of helium with respect 
to hydrogen produces a decrease of the mean molecular 
weight ^ in the envelope as, assuming here a complete 
ionisation : 



M- 



2X 



-Y 



and d^, 



y(_dx + — ). 



In the outer layers of a model computed with microscopic 
diffusion dX > Q, dY < Q and dZ < Q then d^ < as 
Z < X, see Fig. 0. 

This variation of the mean molecular weight produces 
an increase of the radius and consequently a decrease of 
the effective temperature. But, the variation of the metal 
content introduces a concommitant variation of the opaci- 
ty in the radiative regions, which tend to counterbalance 
this effect. 

With the initial abundances given in Table [I], we have 
computed respectively sets of models, i) without diffusion, 
ii) with diffusion of helium only, iii) with diffusion of all 
species. 

The main characteristics of four models of masses 
0.85Mq, 0.8Mq and O.75M0, at the age of lOGyr com- 
puted without and with microscopic diffusion of all species 
are displayed in Table and in Table for the models of 
O.7M0 and O.GMqQ 

During the main-sequence, models with diffusion are 
cooler than models without diffusion; this is essentially 
due to the changes of the mean molecular weight. As 
already stressed by Castellani et al. ( 1997D comparison 

^ There, for the calculations of [Fe/H] with Eq. (m), the differ- 
ential changes between Zfo and Z, due to gravitational settling 
and radiative forces, which act on opposite way, are neglected 
and log(Zpc/Z) = —1.371 is used. 
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Fig. 3. The main-sequence HR evolutionary tracks for 
0.8Mq, [Fe/H] = -0.94, without diffusion (thin full), with 
diffusion of helium only (dash-dot-dash), with diffusion of 
all elements (dashed), and finally with the calibration of 
surface [Fe/H] (heavy full) (the initial metallicity has been 
computed in order to achieve, at 10 Gyr, the initial surface 
metallicity of the other three models). On each evolution- 
ary track the open circles correspond to ages of 3 Gy, 6 Gy 
and 9Gy respectively. For 9Gy, AB is the diffusion shift 
DS, BC the calibration shift CS and AC the global shift 
GS (see text). 



with unphysical models, including diffusion of helium only, 
shows that a significant contribution to the displacement 
in the HR diagram is due to helium settling, in particular 
for low metal content as illustrated on Fig. g. 

In the following, we will call "diffusion shift" DS, the 
translation in the HR diagram from a "standard model" 
to a model including microscopic diffusion, at the same 
mass and the sam,e age. It consists in : 

— a reduction of the effective temperature, decreasing 
with mass from 100 K at 0.85Mq, to almost zero at 
0.6Mq; the stronger increase when [Fe/H] decreases is 
explained by the larger influence of helium stratifica- 
tion at low metallicities; 

— an augmentation of the luminosity of the order of 
0.01 to 0.02 dex, almost independent of [Fe/H], very 
sightly decreasing with mass, almost certainly due to 
the higher temperature and higher helium content of 
the central core of models with diffusion. 

The time dependence of the process shows up through 
the increasing separation between the tracks as evolution 
proceeds. 

5. "Calibration" of the subdvirarfs main-sequence 

Due to the decrease of the surface value of [Fe/H] along 
the evolution, the stars presently observed with a given 
[Fe/H]^ have started their evolution with a larger [Fe/H]Q 
value of their metal content. 
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Table 2. Characteristic of models of O.SSM© (Q), O.SOM© (N), 0.75Mq (R), OJOMq (P) and O.6OM0 (S) at lOGy. 
The models evolved with microscopic diffusion are specified by the labels 'dc' and 'd' respectively for calibrated and 
non-calibrated (see text). Mboi is the bolometric magnitude, (MQboi = 4.75); X^ and Y^, are respectively the surface 
abundances, in mass unit, of hydrogen and helium; Xq and Zq are respectively the initial surface abundances, in mass 
unit, of hydrogen and heavy elements of calibrated models "dc"; [Fc/HJq and [Fe/H]^ are respectively the surface 
values of [Fe/H] at zero age main-sequence and at lOGy; GSt in Kelvin, is the component on the Tes axis of the global 
shift GS (see text) for 10 Gy. 





Q07 


Q09 


Q12 


Q17 


N07 


N09 


N12 


N17 


R07 


R09 


R12 


R17 


Tcs 


5832 


6097 


6401 


6720 


5545 


5822 


6144 


6430 


5222 


5500 


5832 


6116 


Mboi 


4.972 


4.649 


4.243 


3.811 


5.447 


5.179 


4.869 


4.573 


5.903 


5.664 


5.398 


5.163 




Q07d 


Q09d 


Q12d 


Q17d 


N07d 


N09d 


N12d 


N17d 


R07d 


R09d 


R12d 


R17d 


TcS 


5809 


6060 


6322 


6584 


5531 


5801 


6101 


6363 


5215 


5488 


5810 


6080 


Mboi 


4.925 


4.592 


4.174 


3.702 


5.408 


5.136 


4.818 


4.510 


5.869 


5.628 


5.358 


5.119 


X, 


0.802 


0.821 


0.853 


0.880 


0.795 


0.807 


0.825 


0.859 


0.789 


0.797 


0.808 


0.825 


Y. 


0.192 


0.177 


0.146 


0.120 


0.200 


0.190 


0.174 


0.141 


0.206 


0.200 


0.190 


0.175 


[Fe/H]^ 


-0.85 


-1.07 


-1.45 


-2.00 


-0.83 


-1.04 


-1.38 


-1.94 


-0.82 


-1.02 


-1.35 


-1.86 




Q07dc 


Q09dc 


Q12dc 


Q17dc 


N07dc 


N09dc 


N12dc 


N17dc 


R07dc 


R09dc 


R12dc 


R17dc 


Xo 


0.753 


0.755 


0.757 


0.759 


0.753 


0.756 


0.758 


0.759 


0.753 


0.756 


0.757 


0.759 


Zo 


0.007 


0.005 


0.003 


0.001 


0.007 


0.005 


0.002 


0.001 


0.007 


0.004 


0.002 


0.001 


TcS 


5674 


5920 


6202 


6477 


5415 


5672 


5999 


6293 


5117 


5405 


5751 


6054 


Mboi 


5.077 


4.779 


4.382 


3.892 


5.519 


5.262 


4.937 


4.590 


5.955 


5.698 


5.412 


5.144 


X. 


0.797 


0.811 


0.835 


0.872 


0.791 


0.801 


0.816 


0.849 


0.786 


0.795 


0.805 


0.822 


Y* 


0.197 


0.185 


0.163 


0.128 


0.203 


0.195 


0.182 


0.151 


0.208 


0.201 


0.193 


0.117 


[Fe/H]„ 


-0.66 


-0.82 


-1.08 


-1.50 


-0.67 


-0.84 


-1.12 


-1.56 


-0.68 


-0.88 


-1.17 


-1.65 


[Fe/H], 


-0.75 


-0.94 


-1.25 


-1.76 


-0.75 


-0.94 


-1.24 


-1.76 


-0.75 


-0.96 


-1.27 


-1.78 


GSt 


-158 


-177 


-199 


-243 


-130 


-150 


-145 


-137 


-105 


-95 


-81 


-62 



To represent observed data, models have then to 
be "calibrated" to reach the observed surface [Fe/H], 
value at their present age. Assuming a linear dependence 
A [Fe/H] /[Fe/H] with respect to mass and age, for a given 
mass and a given age, one can estimate the initial [Fe/H]Q 
needed to reach the observed value, and start again a new 
evolutionary sequence including diffusion. If needed, this 
process can be iterated. 

This effect adds a new shift of the HR position of the 
models towards the red, called the "calibration shift" CS, 
see Fig. 0. It also increases with age, as it is due to a 
reduction of the surface abundances as evolution proceeds. 
R consists in : 

— a decrease of the effective temperature due to the 
higher initial metal content, which is higher for higher 



metal content and decreases with mass from 140 K at 
0.85Mo, to 30 K at O.GMq, for [Fe/H] = -0.94. 
— a reduction of the luminosity, strongly depending on 
mass and becoming very small for 0.6Mq, which is also 
due to the higher metal content. 

The extend of CS depends on the relative variation of 
[Fe/H] during the evolution, as seen from Table and 
where are displayed the main characteristics of the mod- 
els with enhanced initial [Fe/H]^ (label "dc") in order to 
obtain at the end of the evolution, the surface values of 
[Fe/H], close to the values of Table ||. 

6. The global shift 

We will call "global shift" GS, the combined effect of dif- 
fusion and calibration. Its component in effective tempera- 
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Table 3. Same as Table | for models of 0.70Mq (P) and 0.60Mq (S). 

P07 P09 P12 P17 S07 S09 S12 S17 



Teff 4885 5146 5473 5774 4275 4453 4699 4983 

Mboi 6.354 6.137 5.895 5.688 7.264 7.087 6.883 6.701 



P07d P09d P12d P17d S07d S09d S12d S17d 



Mboi 
[Fe/H] 



4880 5139 

6.325 6.105 

0.784 0.791 

0.211 0.206 



5461 
5.861 
0.798 
0.200 



5756 
5.653 
0.807 
0.193 



4277 
7.239 
0.776 
0.218 



4453 
7.061 
0.782 
0.215 



4699 4983 

6.856 6.673 

0.786 0.789 

0.212 0.211 
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Fig. 4. The main-sequence evolutionary tracks of 0.7 Mq 
models, evolved without diffusion (dashed) and cali- 
brated (full). The solar bolometric magnitude is taken as 

M0boi = 4.75. 



ture is given in Table || and Q for the 10 Gy models. Fig. ^ 
shows, as a function of mass, the variations of effective 




3.6 



Fig. 5. Schematic isochrones at 10 Gy for models with 
[Fe/H]^ = -0.75, -0.94, -1.24 and -1.72; dashed : 
model without diffusion, full : calibrated models with dif- 
fusion. • : locations of stars with observed [Fe/H] within 
— 1.0 < [Fe/H] < —0.45 as derived from Hipparcos data. 
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Fig. 6. Total variation of effective temperature at 10 Gy 
between standard and calibrated models with respect to 
the mass. At 0.85 Mq evolutionary effects become signifi- 
cant. 
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Fig. 7. Mass-luminosity diagram at 10 Gy for standard 
(thin) and calibrated (thick) models. 




Log 10 L/Lq 

Fig. 8. Temperature difference between two isochrones at 
the same luminosity as a function of luminosity, respec- 
tively for 10 Gy (full) and 8Gy (dahed), with [Fe/H]^ = 
— 1.72 - wigles are due to the interpolation procedure. 



tionary track, see Fig. o, due to the time dependence of 
the process. Isochrones are then slightly more vertical. 

Stars with observed [Fe/H] within 

— 1.0 < [Fe/H] < —0.45 and the most accurate fun- 
damental parameters 



studied by Lebreton et al. (1997) 
have been added to Fig. ^. As already said, though micro- 
scopic diffusion tends to reduce the discrepancy between 
the effective temperature of observed and theoretical 
models, the effect is not sufficient to solve the difficulty. 

For [Fe/H] = —1.72, Fig. || shows, as a function of lu- 
minosity and effective temperature, the effective temper- 
ature difference between two isochrones at the same lumi- 
nosity : the standard one experiences no diffusion and its 
initial abundances are equal to the observed ones, while 
the second one is computed using our proposed calibra- 
tion. 



temperature for the three different metal content consid- 
ered here. 

Both diffusion and calibration reduce the effective tem- 
perature, but in an opposite way with respect to [Fe/H]. 
The effect of diffusion increases while the effect of calibra- 
tion decreases with [Fe/H]. 

In luminosity the situation is different. The calibration 
procedure compensates partly the luminosity increase due 
to diffusion. The variation can reach 0.03 dex, depending 
both on mass and [Fe/H]. The mass dependence remains 
strong, from 0.06 dex at 0.85Mq, to 0.01 dex at 0.6Mq. 
It produces a distortion in the Mass-Luminosity relation, 
see Fig. 0, depending on [Fe/H], an effect which could be 
detected with high quality data on a reasonable set of 
objects. 

The displacement in the HR diagram is seen as a trans- 
lation towards lower effective temperatures, see Fig. ^. It 
produces a slight variation of the curvature of the evolu- 



7. Discussion 

We have estimated in this paper the influence of gravita- 
tional and thermal settling, on the global parameters of 
low main-sequence stars. As we have introduced only mi- 
croscopic flows in the computations except in convective 
zones, the stratification is probably overestimated. 

We have shown that this physical process, if at work 
permanently, produces an effective temperature shift to- 
wards the red. For population II stars it exceeds generally 
100 K. 

Compared to the most reliable data available at 
present, the predicted shift is of the right order of magni- 
tude, but probably slighly unsufficient to eliminate com- 
pletely the discrepancy shown by Lebreton et al. (1997). 



However, as already mentioned, other effects could also be 
responsible for redshifts of that kind. So that more precise 
characteristic behaviors are needed to decipher among the 
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different possibilities. Let us cite a few signatures of the 
microscopic diffusion. 

— the predicted shift increases with age and decreases 
when [Fe/H] increases; 

— the mass-himinosity relationship is modified depending 
on Z as at a given mass and age, as diffused and cali- 
brated models are less luminous than standard ones; 

— the slope of isochrones is slightly increased. 

Though very important improvement on the quality of 
the observational data on the nearest members of the old 
population have been obtained recently, they remain too 
coarse and too few to allow unambiguous tests. More ob- 
jects with accurate fundamental parameters are needed to 
validate definitely this hypothesis. 

A lot of prog ress is expected from the GAIA mission 
(Ferryman et al. 1997b ), which will provide even more ac- 
curate distances and sometimes masses for a larger num- 
ber of stars. But improvements in the determination of 
atmospheric parameters, as well as in methods of com- 
paring observed and predicted quantities are also badly 
needed. 
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